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Abstract 

Primary and atmospheric cosmic-ray spectra were precisely measured with the 
BESS-TeV spectrometer. The spectrometer was upgraded from BESS-98 to achieve 
seven times higher resolution in momentum measurement. We report absolute fluxes 
of primary protons and helium nuclei in the energy ranges, 1-540 GeV and 1- 
250 GeV/n, respectively, and absolute flux of atmospheric muons in the momentum 
range 0.6-400 GeV/c. 
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1 Introduction 



The absolute flux and spectral shape of primary cosmic rays are the basis to 
discuss the origin and the propagation history of the cosmic rays in the Galaxy. 
The spectrum is also essential as an input to calculate spectra of cosmic- 
ray antiprotons and positrons which are secondary products of cosmic-ray 
interactions with the interstellar gas. Recently, the importance of the cosmic- 
ray spectrum has been emphasized in connection with the study of neutrino 
oscillation observed in atmospheric neutrinos [1]. In contrast with a long- 
baseline experiment, such as K2K [2], the feature of the oscillation study with 
atmospheric neutrinos lies in its wide energy range. For example, the Super 
Kamiokande water Cherenkov detector [3] can observe neutrinos from 0.1 to 
100 GeV. In order to estimate an accurate flux of atmospheric neutrinos up 
to around 100 GeV, we need to know the primary cosmic-ray flux well above 
100 GeV and the feature of hadronic interactions in that energy region. It is 
crucial to measure primary cosmic-ray flux by determining the absolute energy 
of particles. Measurements of atmospheric muon spectrum is also important 
to check and improve our understanding of hadronic interactions. 

The Balloon-borne Experiment with a Superconducting Spectrometer (BESS) 
[4,5] has been carried out since 1993 to perform highly sensitive searches for 
cosmic-ray antiparticles and precise measurements of absolute flux of various 
cosmic-ray components. The absolute fluxes of primary protons and helium nu- 
clei were precisely measured in the energy ranges, 1-120 GeV and 1-54 GcV/n, 
respectively, by a balloon observation in 1998 [6]. The overall uncertainties in 
the measurements were less than 5 % for protons and 10 % for helium nuclei. 
The absolute flux of atmospheric muons was measured from 0.6 to 20 GeV/c 
at sea level [7] and from 0.6 to 100 GeV/c at a mountain altitude [8]. The 
overall uncertainty in the measurements was less than 10 %. In order to ex- 
tend the energy range of the precise measurements of cosmic-ray flux up to 
higher energy, the BESS-TeV spectrometer was developed. A magnetic rigid- 
ity (i? = pc/Ze) of an incident particle was accurately measured by a new 
tracking system. The maximum detectable rigidity (MDR) was signiflcantly 
improved from 200 GV to 1.4 TV. 

Cosmic-ray observations were carried out with the BESS-TeV spectrometer at 
a balloon altitude and at sea level in 2002. We report precise measurements of 
absolute fluxes of primary protons and helium nuclei, and atmospheric muons. 
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As to primary cosmic rays, we focused on a energy range above 1 GeV/n in this 
paper. The BESS spectrometer has measured primary cosmic-ray flux down 
to 0.2 GeV/n, which provides significant information on the solar modulation 
effect and the propagation of cosmic rays [9] . Measurements of the low-energy 
cosmic-ray flux with the BESS-TeV spectrometer will be reported elsewhere 
together with a series of BESS balloon- flight data [10]. 



2 The BESS-TeV spectrometer 



The BESS spectrometer [11] was upgraded to achieve a significantly high rigid- 
ity resolution. The upgraded spectrometer, "BESS-TeV" [12], was equipped 
with newly-developed drift chambers. 

As shown in Fig. 1, the detector has a unique feature of a cylindrical con- 
figuration realized by a thin superconducting solenoid [13]. The configuration 
resulted in a large and almost constant geometrical acceptance and a uniform 
detector performance for various incident angles and positions. 

In the central region, the solenoid with a diameter of 1 m provides a uniform 
magnetic field of 1 T. The field variation is less than 2.5 % along a typi- 
cal trajectory of an incoming particle. A defiection (R^^) of the trajectory is 
measured by a central jet-type drift chamber (JET), two inner drift cham- 
bers (IDCs) and two outer drift chambers (ODCs), all of which were newly 
developed for the BESS-TcV spectrometer. Inside JET and IDCs a trajectory 
was determined by simple circular fitting [14] using up to 52 hit points. Each 
hit point was measured with a spacial resolution of 150 fim. More accurate 
defiection was measured by adding 8 hit points inside ODCs which are placed 
outside the solenoid at a radius of 0.8 m. The hit positions inside the ODCs 
were used above 100 GV, where the effect of multiple scattering in the detec- 
tor material is negligibly small. Fig. 2 shows a distribution of the defiection 
resolution (Ai?^^) obtained with all the drift chambers. The defiection res- 
olution was evaluated in the track-fitting procedure for cosmic-ray protons. 
Those of other spectrometers used in previous balloon experiments [6,15,16] 
are also shown. Each area of the histogram is normalized to unity. The peak 
position of 0.7 TV~^ corresponds to the MDR of 1.4 TV. In order to pro- 
vide an absolute reference position for their calibration, a scintillating fiber 
counter system (SciFi) was attached to the top and the bottom walls of the 
ODCs. SciFi consists of two layers of 1 mm x 1 mm square-shaped scintil- 
lating fibers and covers the central cell of each ODC. The scintillating fiber 
layers are shifted by 0.5 mm with each other. SciFi can measure a hit position 
with an accuracy of the overlap width of two fiber layers, i.e. 0.5 mm. With a 
sufficiently large number of events, the track position can be determined with 
a better spatial resolution than that of ODCs. 
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Time-of-flight (TOF) hodoscopes [17] provide the velocity and energy 
loss (dii^/d,x) measurements. A j3^^ resolution of 1.4 % was achieved in the 
experiment. The data acquisition sequence is initiated by a first-level TOF 
trigger, which is a simple coincidence of signals from the upper and lower 
TOF counters. The trigger efficiency was evaluated to be 99.4±0.2 % by a 
secondary proton beam at KEK 12 GeV proton synchrotron. The trigger rates 
were 1 kHz and 30 Hz at a balloon altitude and sea level, respectively. During 
the balloon observation, one out of every ten events were recorded to sample 
unbiased trigger events. An auxiliary trigger is generated by a signal from 
a Cherenkov counter [18] to record particles above threshold energy without 
bias or sampling. The efficiency of the Cherenkov trigger was evaluated as the 
ratio of the number of Cherenkov-triggered events to the unbiased triggered 
events. It was found to be 94. 1± 2.0% and 97.3±4.0 % for relativistic protons 
and helium nuclei, respectively. For the determination of primary proton and 
hehum fluxes, the Cherenkov-triggered events were used above 10 GeV and 
5.5 GeV/n, respectively, and the TOF-triggered events were used below these 
energies. During the ground observation, all the TOF-triggered events were 
recorded to determine the atmospheric muon flux. 



3 Observations 



The BESS-TeV spectrometer was launched by a balloon from Lynn Lake, 
Manitoba, Canada (56.5°N, 101. 0°W), on 7th August 2002. After about four- 
hour ascending, the payload reached a floating altitude of 37 km (residual 
atmosphere of 4.8 g/cm^). The geomagnetic cutoff rigidity was 0.5 GV or 
smaller throughout the flight. The total live time of the data-taking was 38,215 
seconds (10.6 hours) during the floating period. Among them 11 % (1.2 hours) 
of data taken around the sun rise were not used in this analysis, because the 
rapid temperature variation might introduce a large systematic error in the 
chamber calibration. 

The ground observation was carried out at KEK located at Tsukuba, Japan 
(36.2°N, 140.1°W), during a period of lst-6th October 2002. Tsukuba is lo- 
cated 30 m above sea level with the vertical cutoff rigidity of 11.4 GV. The 
variations in the atmospheric pressure and temperature during the observa- 
tion are shown in Fig. 3. The atmospheric temperature data were obtained 
from Ref. [19]. The arrows indicate the data-taking periods used to determine 
the atmospheric muon flux. We did not use the data during a period where 
the variation in atmospheric pressure was large. The mean (root-mean-square) 
atmospheric pressure and temperature during the periods used for the flux de- 
termination were 1032.2 g/cm^ (4.4 g/cm^) and 20.8 °C (3.7 °C), respectively. 
The total live time of the data-taking was 329,403 seconds (91.5 hours). 
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4 Data analysis 



4.1 Event reconstruction 



Each hit position inside the drift chambers was calculated from the drift time 
digitized by a flash analog-to-digital converter. The calculation was carried 
out based on a relation between the hit position and the drift time {x-t re- 
lation). The x-t relation was precisely calculated by a drift chamber sim- 
ulation package, GARFIELD [20], and a gas property simulation package, 
MAGBOLTZ [21]. Although the chambers were constructed carefully with a 
tolerance of 100 /xm, there was a small position deviation of wires and fleld- 
shaping patterns, which could locally modify the electric field. In order to take 
account of the limited accuracy in the chamber manufacturing, a correction 
was commonly applied to the calculated x-t relation throughout the experi- 
ments. The correction was obtained to minimize the in the fitting of straight 
tracks of clean muon events observed on the ground without magnetic field. 
The correction was as small as expected from the accuracy in the chamber 
manufacturing. During the observations, the x-t relation was affected by the 
variation in the pressure and temperature of the chamber gas. In order to take 
account of these time-dependent variations, the x-t relation was calibrated for 
each data-taking run. Especially in calibrating the x-t relation of ODCs, an 
absolute reference positions were provided by SciFi, which are not affected by 
the variation in the pressure nor temperature. 

The same calibration procedure was applied to the muon events observed on 
the ground without magnetic field, where the straight tracks were used as an 
absolute reference of events with infinite rigidity. We checked the defiection of 
the muon events by applying circular-fitting to the calibrated hit points inside 
the JET and IDCs. The mean value of the obtained deflection was smaller 
than (10 TV)~^. Therefore, the systematic shift in the deflection measurement 
originated from this calibration procedure should be smaller than (10 TV)~^. 

The precise alignment of the ODCs with respect to the JET was determined by 
checking the consistency between tracks reconstructed by JET and hit points 
measured by ODGs. In order to minimize the effect of multiple scattering in 
the detector material, we selected events whose rigidity was measured to be 
higher than 4 GV by JET and IDCs. The chamber alignment was calibrated 
for each run. We checked the variation in the calibrated position of ODCs in 
84 runs during the ground observation. The root-mean-square of the variation 
was smaller than 20 /xm. This variation might introduce a systematic error in 
the deffection measurement of AR~^ — (5.3 TV)~^. 
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4-2 Event selection 



We selected events with a single track fully contained inside the fiducial volume 
defined by the central six columns out of eight columns in the JET. It ensured 
that the track is long enough for reliable rigidity measurement. In order to 
obtain a nearly vertical flux of atmospheric muons, an additional cut on the 
zenith angle {9z) was applied as cos^^ > 0.98 below 20 GV and cos^^ > 0.90 
above 20 GV. We have checked the change of the atmospheric muon flux above 
20 GV by requiring two conditions, cos^z > 0.98 and cos^^^ > 0.90. We found 
that the requirement to zenith angle could be relaxed to cos6'z > 0.90 without 
changing the observed flux. For protons and helium nuclei zenith angle cut was 
not applied. Due to the detector geometry, however, zenith angle is limited as 
cos^z > 0.80 above 2 GV. 

A single-track event was deflned as an event which has only one isolated track 
inside the JET and one or two hit counters in each layer of the TOF ho- 
doscopes. The single-track selection eliminated rare interacting events. To es- 
timate the efficiency of the single-track selection, Monte Carlo simulations with 
GEANT3/4 [22,23] were performed. The probabihty that each particle could 
pass through the selection was evaluated by applying the same selection crite- 
ria to the Monte Carlo events. The resultant efficiency of the single-track event 
selection for protons was 82.0±2.6 % at 20 GV and 79.9±4.0 % at 200 GV, 
the efficiency for helium nuclei was 69.6±6.6 % at 20 GV and 66.9±7.2 % at 
200 GV, and the efficiency for muons was 96.8±1.5 % at 20 GV and 95.9±1.5 % 
at 200 GV. In order to ensure the track reconstruction in the z-direction (per- 
pendicular to the bending plane), we required a consistency of hit position 
in the 2;-direction measured by two independent detectors; the drift chamber 
and the TOF hodoscope. The efficiency of this selection was estimated to be 
97.2±1.0 % for protons, 95.5±1.5 % for helium nuclei, and 98.1±1.0 % for 
muons. 

Particle identification was performed by requiring proper dE/dx measure- 
ments with both upper and lower layers of the TOF hodoscopes and as 
functions of rigidity. Figs. 4 and 5 show the selection criteria for protons. He- 
lium nuclei and muons were identified in the same manner. The efficiencies of 
dE/dx selection were estimated with another sample selected by independent 
measurement of energy loss inside the JET. We found that 96.3±0.3 % of 
protons, 92.3±0.9 % of helium nuclei and 99.9±0.1 % of muons were properly 
identified. Since the f3~^ distribution is well described by Gaussian and a half- 
width of the (3^'^ selection band was set at 4 a, the efficiency is very close to 
unity. Following this particle identification procedure, 444,578 proton, 38,006 
helium, and 688,983 muon candidates were obtained. 



6 



4-3 Backgrounds 



As shown in Fig. 5, protons were identified without contamination by exam- 
ining f3~^ distribution below 2 GV. Above 2 GV, there were contamination of 
hght particles such as positrons, pions and muons. The ratio of hght particles 
to protons was observed to be 4.5 % at 1 GV and 1.1 % at 2 GV, which was 
expressed by a power law with an index of —2.1. This tendency was repro- 
duced well by a Monte Carlo simulation [24] based on the DPMJET-111 event 
generator [25]. Above 10 GV, the simulated n'^/p ratio shows almost con- 
stant value of 0.2 %. In this analysis the contamination of the light particles 
were subtracted by assuming the power law, which would underestimate the 
background of light particles. However, The contamination is smaller than the 
statistical errors. Above 3 GV, there were contamination of deuterons, which 
was observed to be 2 % at 3 GV. According to our previous measurement [9], 
the d/p flux ratio decreases with increasing rigidity. The ratio should follow 
a decrease in escape path lengths of primary cosmic-ray nuclei [26]. Since the 
contamination of deuterons was as small as the statistical errors, and should 
decrease with increasing rigidity, no subtraction was made for the deuteron 
contamination. Therefore, above 3 GV hydrogen nuclei were selected, which 
included a small amount of deuterons. Helium nuclei were identified clearly by 
redundant charge measurement inside the upper and lower TOF hodoscopes. 
The helium candidates contained both ^He and ^He. In conformity with pre- 
vious experiments, all doubly charged particles were treated as ^He. 



Among the muon candidates observed on the ground, there were contamina- 
tions of electrons, positrons and protons. In our previous work, the ratio of 
electrons and positrons to muons was measured to be 1.5 % at 0.5 GV and 
0.3 % at 1 GV with an electro-magnetic shower counter [7]. Since the BESS- 
TeV spectrometer is not equipped with an electro- magnetic shower counter, 
no subtraction was made for the contamination of electrons and positrons. 
However, the contamination was smaller than the statistical errors. Protons 
were identified by examining distribution below 2.5 GV. The p//x fiux 
ratio was observed to be 3.3 % at 1 GV and 1.2 % at 2.5 GV, which was 
expressed by a power law with an index of —1.0. Below 7 GV, the power law 
agreed well with a previous measurement of the p/fi ratio [27], which shows, 
however, constant value of about 0.4 % above 7 GV. The proton background 
was subtracted from muon candidates assuming the power law below 7 GV 
and constant value of 0.4 % above 7 GV. 
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4-4 Normalization and corrections 



In order to obtain the absolute flux of protons, helium nuclei and muons 
at the top of the detector, energy loss by ionization inside the detector ma- 
terial, live time and geometrical acceptance were estimated. The energy of 
each incoming particle was calculated by integrating the energy losses inside 
the detector tracing back along the particle trajectory. The total live time of 
the data-taking was measured exactly by counting 1 MHz clock pulses with 
a scaler system gated by a "ready" status that controls the first-level trig- 
ger. The geometrical acceptance defined for this analysis was calculated as 
a function of rigidity with a simulation technique [28]. The geometrical ac- 
ceptance is 0.0886±0.0003 m^sr for protons and helium nuclei at 10 GV and 
0.0302±0.0001 m^sr for muons at 10 GV. The acceptance for muons is about 
1/3 of that for protons and helium nuclei because of the additional require- 
ment on the zenith angle (cos^z > 0.98) as described in Section 4.2. The simple 
cylindrical shape and the uniform magnetic field make it simple and reliable to 
determine the geometrical acceptance precisely. The error on the acceptance 
calculation was estimated to be 0.3 % which represents the uncertainty of the 
detector alignment 

In order to obtain the absolute flux of primary protons and helium nuclei at the 
top of the atmosphere, interaction loss and secondary particle production in 
the residual atmosphere were estimated. According to the Monte Carlo studies, 
probabilities for primary cosmic rays to penetrate the residual atmosphere of 
4.8 g/cm^ are 93.8±0.7 % and 91.3±2.0 % for protons and helium nuclei, 
respectively, at 10 GeV and almost constant over the entire rigidity range 
discussed here. Atmospheric secondary protons were subtracted based on a 
calculation for the maximum solar activity epoch by Papini et al. [29]. A 
secondary-to- primary proton ratio is 4.8±0.5 % at 1 GeV and 1.7±0.2 % above 
10 GeV. Atmospheric secondary helium nuclei above 1 GeV/n are dominated 
by the fragments of heavier cosmic-ray nuclei (mainly carbon and oxygen). 
The flux ratio of the atmospheric secondary helium to the primary carbon and 
oxygen was estimated to be 14 % at a depth of 4.8 g/cm^, based on the total 
inelastic cross sections of CNO-|-Air interactions and the hehum multiphcity in 
^^C -|- CNO interactions [30]. The total correction of atmospheric secondary 
helium produced by all the primary nuclei with Z>2 was estimated to be 
1.6±0.5 % at 1 GeV/n and 2.1±0.6 % at 10 GeV/n. 



4-5 Spectrum deformation effect 

Because of the limited accuracy of the rigidity measurement and the steep 
spectral shape, the observed spectrum may suffer deformation. The limited 
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accuracy of the rigidity measurement was decomposed into two sources: (i) a 
finite resolution in the rigidity measurement, and (ii) a small shift in the 
measured deflection which is caused by the misalignment of the ODCs. We 
estimated the systematic error in the spectrum measurement caused by such 
limited accuracy of the rigidity measurement. The systematic error was esti- 
mated by a Monte Carlo simulation. In the simulation, detailed response of 
the drift chamber was implemented, such as the distribution of primary ion- 
ization clusters, the diffusion in the gas, the fluctuation of avalanche gain, and 
the digitization in the readout electronics. The simulated rigidity resolution 
reproduced well the experimental result shown in Fig. 2. Therefore the effect 
of spectrum deformation should be correctly estimated by this simulation. As 
an input spectrum, a power-law function was used with a spectral index of 
—2.7 for both primary protons and helium nuclei. For atmospheric muons, the 
spectral index was chosen as —3.2, which was obtained by fitting a power-law 
function to the muon spectrum obtained with this experiment in the rigid- 
ity range 100-400 GV. Fig. 6 shows the effect of spectrum deformation as a 
ratio of simulated spectrum to the input spectrum. The open squares show 
the spectrum deformation when ODCs are aligned correctly. In this case, the 
deformation is caused only by the effect of (i). The change of the spectrum 
is less than 5 % for protons and helium nuclei below 1 TV, and for muons 
below 400 GV. The upward and downward triangles show the spectrum de- 
formation when ODCs arc artificially displaced by ±20;um. We found that the 
shift of the measured dcficction was smaller than (5.3 TV)^^ as described in 
Section 4.1. The change of the spectrum was less than 2 % below 100 GV and 
10 % at 500 GV for primary protons and helium nuclei, and less than 3 % 
below 100 GV and 13 % at 400 GV for nmons. Since the input spectrum of 
atmospheric muons is steeper than that of primary protons and helium nuclei, 
the spectrum deformation effect is more significant for atmospheric muons. 
We did not apply any corrections, such as unfolding, on the observed spectra. 
The spectrum deformation effect estimated above was treated as a systematic 
error on the observed spectra, which was as small as the statistical error over 
the entire rigidity range in this analysis. 



5 Results and discussions 



We have obtained the absolute fiuxes of primary protons in the range 1- 
540 GeV and hehum nuclei in the range 1-250 GeV/n at the top of the atmo- 
sphere from the BESS-TeV balloon- flight data in 2002. We have obtained the 
absolute flux of muons in the range 0.6-400 GcV/c at sea level (30 m a.s.l.) 
from the BESS-TeV ground-observation data in 2002. The results of protons 
and helium nuclei are summarized in Tables 1 and 2, respectively, and the 
results of atmospheric muons are summarized in Table 3. The overall uncer- 
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tainties including both statistical and systematic errors were less than ±15 % 
for protons, ±20 % for helium nuclei, and ±20 % for muons. 

The results of primary proton and helium spectra are shown in Fig. 7 in com- 
parison with other experiments with magnetic spectrometers [6,15,16,31,32,33]. 
Discrepancies in the observed spectra below 10 GeV for protons and 5 GeV/n 
for hehum nuclei come from the difference of solar activity (around minimum 
in 1998 and around maximum in 2002). The detailed analysis of the solar 
modulation affecting the low-energy spectra in a series of BESS flights will be 
discussed elsewhere [10]. 

Although a small systematic shift of 2-3 % was found in absolute proton 
flux between the results of BESS-TeV and BESS-98 from 30 to 100 GeV, the 
results are well consistent within the overall uncertainty of 5 %. The main 
source of this uncertainty was the aerogel trigger efficiency. The accuracy of 
the efficiency estimation was limited by the statistics of the event sample to 
be 2.0 % and 3.0 % for BESS-TeV and BESS-98, respectively. 

Our resultant spectral shape of protons and helium nuclei is very similar to 
that measured by AMS [31,32]. Above 30 GeV, the absolute proton flux mea- 
sured by BESS-TeV and AMS shows good agreement within 5 %. However, 
there is 15 % discrepancy in the absolute flux of helium nuclei. Both pro- 
ton and helium spectra by CAPRICE-98 [15] show steeper spectra than our 
results. 

At high energies, the spectrum F may be parameterized by a power law in 
kinetic energy, E]^, as F — ^E^'^. The fitting range was chosen to be 30- 
540 GeV for protons and 20-250 GeV/n for helium nuclei so that the solar 
modulation effect was negligible. The best fit values and uncertainties for 
protons ($p and 7p) and helium nuclei ($He and 7He) were obtained as 

% = (1.37 ± 0.06(sta.) ± 0.11(sys.)) x 10^ (m^ • sr ■ s • GeV)"^ 

7p = 2.732 ± 0.011(sta.) ± 0.019(sys.) 

(7.06 ± 0.94(sta.) ± 1.17(sys.)) x 10^ (m^ ■ sr • s • (GeV/n))"^ 
2.699 ± 0.040(sta.) ± 0.044(sys.) 

respectively. The two parameters were strongly correlated, with a correlation 
coefficient of —0.98. 

The result of atmospheric muon spectrum is shown in Fig. 8 in comparison 
with other absolute flux measurements by magnetic spectrometers [7,8,34,35]. 



and 

7He = 
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The discrepancy below 30 GeV/c among the observed muon spectra is mainly 
due to the difference in altitudes. Fig. 9 shows charge ratios of atmospheric 
muons observed in a series of muon measurement with the BESS experiments. 
The difference in the charge ratios observed in Japan and Canada comes from 
the different geomagnetic cutoff rigidity [7]. 



6 Conclusion 

We have measured energy spectra of primary protons in the range 1-540 GeV 
and helium nuclei in the range 1-250 GeV/n by a balloon observation, and 
momentum spectrum of atmospheric muons in the range 0.6-400 GeV/c by a 
ground observation at sea level. 

The overall uncertainties were less than 15 % for protons, less than 20 % for 
helium nuclei, and less then 20 % for muons. Primary cosmic-ray spectra pro- 
vide fundamental information on the origin and propagation history of the 
cosmic rays. The results also provide accurate input spectra to predict an 
atmospheric neutrino flux. The result of atmospheric muon spectrum will im- 
prove the accuracy of atmospheric neutrino calculation in the wide momentum 
range. 
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Deflection resolution AR ^ (TV) 

Fig. 2. Distribution of the deflection resolution for protons evaluated for each event 
in the track-fitting procedure. Each area of the histogram is normalized to unity. 
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Fig. 3. Atmospheric pressure and temperature during the ground observation. Ar- 
rows indicate the data-taking periods used to determine the atmospheric muon flux. 
The atmospheric temperature data were obtained from Ref. [19] . 
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Fig. 4. Proton band in dE/dx (top TOF) vs. rigidity obtained from the balloon 
observation. A dE/dx in the bottom TOF was also checked. The superimposed 
graph shows the selection criteria for protons and helium nuclei above 10 GV. 
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1 10 10^ 10^ 

Rigidity (GV) 

Fig. 5. Scatter plot of j3~^ vs. rigidity obtained from the balloon observation after 
proton dE/dx selection. The superimposed graph shows the selection criteria for 
protons above 10 GV. 



18 



1.1 



el 
o 



0.9 



(a) p,He 



This analysis 



A A 



AAA 



-K-{1I-'[I--II-'B'-II-H'-H-V'-III-III-H' 



Finite resolution only 
Alignment -20 \im 
Alignment +20 \im 




10 



10 10 
Rigidity (GV) 



Fig. 6. Spectrum deformation effects for (a) primary protons and helium nuclei and 
(b) atmospheric muons. The deformation effects are shown as ratios of simulated 

spectra (Fsim) to the input spectra (Fini) as functions of rigidity. Open squares show 
the case when ODCs are aligned correctly. Upward and downward triangles show 
the case when ODCs are artificially displaced by ±20//m. 
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Fig. 7. Absolute differential energy spectra of primary protons and helium nuclei 
multiplied by E^-^. The spectra obtained by other experiments [6,15,16,31,32,33] 
are also shown. 
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site altitule 
BESS-TeV2002 Tsukuba 30 m 

BESS 1999 Mt. Norikura 2770 m 
BESS 1995 Tsukuba 30 m 

BESS 1997/98/99 Lynnn Lake 360 m 
CAPICE 1997 Ft. Sumner 1270 m 
CAPICE1994 LymiLake 360 m 
MASS 1989 Prince Albert 600 m 
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Fig. 8. Absolute differential momentum spectrum of atmospheric muons multiplied 
by P^-^. The spectra obtained by other experiments [7,8,34,35] are also shown. The 
geomagnetic cutoff rigidity at each site is indicated as Rc- The difference below 
30 GeV/c among the spectra is mainly due to the different altitudes. 
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Fig. 9. Charge ratios of atmospheric muons. Only statistical errors are included. The 
charge ratios obtained by other experiments [7,8] are also shown. The geomagnetic 
cutoff rigidity at each site is indicated as Re- The difference among the charge ratios 
is mainly due to the difference in geomagnetic cutoff rigidity. 
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Table 1 

Primary proton flux at the top of the atmosphere. 



Energy range 
(GeV) 


Ek 
(GeV) 


FluxiAFsta ± AFsys 

(m-2sr-is-iGeV-i) 


1.00- 


1.17 


1.08 


3.50 ±0.03 


±0.11 


xlO^ 


1.17- 


1.36 


1.26 


3.22 ±0.03 


±0.10 


xlO^ 


1.36- 


1.58 


1.47 


2.98 ±0.02 


±0.09 


xlO^ 


1.58- 


1.85 


1.71 


2.71 ±0.02 


±0.08 


xlO^ 


1.85- 


2.15 


2.00 


2.41 ±0.02 


±0.07 


xlO^ 


2.15- 


2.51 


2.33 


2.08 ±0.02 


±0.06 


xlO^ 


2.51- 


2.93 


2.71 


1.74 ±0.01 


±0.05 


xlO^ 


2.93- 


3.42 


3.16 


1.45 ±0.01 


±0.04 


xlO^ 


3.42- 


3.98 


3.69 


1.19 ±0.01 


±0.03 


xlO^ 


3.98- 


4.64 


4.30 


9.52 ±0.08 


±0.27 


xlO^ 


4.64- 


5.41 


5.01 


7.35 ±0.07 


±0.21 


xlOi 


5.41- 


6.31 


5.84 


5.63 ±0.05 


±0.16 


xlO^ 


6.31- 


7.36 


6.81 


4.34 ±0.04 


±0.12 


xlO^ 


7.36- 


8.58 


7.93 


3.15 ±0.03 


±0.09 


xlO^ 


8.58- 


10.0 


9.25 


2.25 ±0.03 


±0.06 


xlO^ 


10.0- 


11.7 


10.8 


1.59 ±0.01 


±0.05 


xlOi 


11.7- 


13.6 


12.6 


1.12 ±0.01 


±0.04 


xlOi 


13.6- 


15.8 


14.7 


7.71 ±0.04 


±0.28 




15.8- 


18.5 


17.1 


5.33 ±0.03 


±0.20 




18.5- 


21.5 


19.9 


3.63 ±0.02 


±0.14 




21.5- 


25.1 


23.2 


2.48 ±0.02 


±0.10 




25.1- 


29.3 


27.1 


1.62 ±0.01 


±0.06 




29.3- 


34.1 


31.6 


1.09 ±0.01 


±0.04 




34.1- 


39.8 


36.8 


7.17 ±0.08 


±0.30 


xlO-i 


39.8- 


46.4 


42.9 


4.84 ±0.06 


±0.21 


xlO-i 


46.4- 


54.1 


50.0 


3.15 ±0.05 


±0.14 


xlO-i 


54.1- 


63.1 


58.3 


2.07 ±0.03 


±0.09 


xlO-i 


63.1- 


73.6 


68.0 


1.34 ±0.03 


±0.06 


xlO-i 


73.6- 


85.8 


79.2 


9.09 ±0.19 


±0.43 


xlO-2 


85.8- 


100. 


92.3 


5.75 ±0.14 


±0.28 


xlQ-^ 


100.- 


126. 


112. 


3.43 ±0.11 


±0.18 


xlQ-^ 


126.- 


158. 


140. 


1.98 ±0.07 


±0.11 


xlO-2 


158.- 


200. 


177. 


1.00 ±0.05 


±0.06 


xlO-2 


200.- 


251. 


222. 


5.42 ±0.31 


±0.34 


xlO-3 


251.- 


316. 


281. 


2.46 ±0.19 


±0.17 


xlO-3 


316.- 


398. 


352. 


1.62 ±0.14 


±0.14 


xlO-3 


398.- 


541. 


463. 


7.47 ±0.69 


±0.78 


xlO-4 
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Table 2 

Primary helium flux at the top of the atmosphere. 



Energy range FluxiAFgta ± AF^ 



(GeV/n) 


(GeV/n) 


(m 


2 —1 
sr 


-i(GeV/n)- 


) 


1.00- 


1.17 


1.08 


5.22 


±0.14 


±0.36 


xlOi 




1.17- 


1.36 


1.26 


4.78 


±0.13 


±0.33 


xlOi 




1.36- 


1.58 


1.47 


4.02 


±0.11 


±0.28 


xlQi 




1.58- 


1.85 


1.71 


3.21 


±0.09 


±0.22 


xlQi 




1.85- 


2.15 


2.00 


2.62 


±0.08 


±0.18 


xlO^ 




2.15- 


2.51 


2.33 


2.17 


±0.06 


±0.15 


xlOi 




2.51- 


2.93 


2.71 


1.81 


±0.05 


±0.12 


xlOi 




2.93- 


3.42 


3.16 


1.37 


±0.04 


±0.09 


xlOi 




3.42- 


3.98 


3.69 


9.77 


±0.35 


±0.67 






3.98- 


4.64 


4.29 


7.67 


±0.28 


±0.53 






4.64- 


5.41 


4.98 


5.71 


±0.23 


±0.39 






5.41- 


6.31 


5.84 


3.98 


±0.06 


±0.32 






6.31- 


7.36 


6.80 


2.83 


±0.04 


±0.22 






7.36- 


8.58 


7.94 


2.07 


±0.03 


±0.16 






8.58- 


10.0 


9.24 


1.48 


±0.03 


±0.12 






10.0- 


11.7 


10.8 


1.02 


±0.02 


±0.08 






11.7- 


13.6 


12.6 


6.76 


±0.16 


±0.54 


xlO^ 


1 


13.6- 


15.8 


14.7 


4.71 


±0.12 


±0.38 


xlO" 


1 


15.8- 


18.5 


17.0 


3.27 


±0.10 


±0.26 


xlO- 


1 


18.5- 


21.5 


19.9 


2.13 


±0.07 


±0.17 


xlO" 


1 


21.5- 


25.1 


23.2 


1.46 


±0.05 


±0.12 


xlO" 


1 


25.1- 


29.3 


27.1 


9.67 


±0.41 


±0.79 


xlO" 


2 


29.3- 


34.1 


31.4 


5.89 


±0.30 


±0.48 


xlO^ 


2 


34.1- 


39.8 


36.7 


4.23 


±0.24 


±0.35 


xlO" 


2 


39.8- 


46.4 


42.9 


2.85 


±0.18 


±0.24 


xlO" 


2 


46.4- 


54.1 


49.9 


1.84 


±0.13 


±0.15 


xlO" 


2 


54.1- 


73.6 


62.5 


9.40 


±0.92 


±0.90 


xlO- 


3 


73.6- 


100. 


86.1 


4.14 


±0.53 


±0.41 


xlO" 


3 


100.- 


136. 


116. 


2.16 


±0.33 


±0.22 


xlO" 


3 


136.- 


251. 


175. 


5.53 


±0.92 


±0.64 


xlO" 


4 
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Table 3 

Atmospheric muon flux at sea level. 



Momentum 

range 



(GcV/c) 


(GeV/c) 


(m 


9 — 1 

^sr -^s 


-i(GeV/c)-: 


0.576- 


0.621 


0.599 


1.34 


±0.02 


±0.03 


xlO^ 


0.621- 


0.669 


0.645 


1.33 


±0.02 


±0.03 


xlO^ 


0.669- 


0.720 


0.694 


1.32 


±0.02 


±0.03 


xlO^ 


0.720- 


0.776 


0.748 


1.27 


±0.02 


±0.03 


xlQi 


0.776- 


0.836 


0.806 


1.26 


±0.02 


±0.03 


xlO'^ 


0.836- 


0.901 


0.869 


1.23 


±0.01 


±0.03 


xlO^ 


0.901- 


0.970 


0.934 


1.21 


±0.01 


±0.03 


xlOi 


0.970- 


1.04 


1.01 


1.16 


±0.01 


±0.02 


xlOi 


1.04- 


1.13 


1.08 


1.10 


±0.01 


±0.02 


xlO^ 


1.13- 


1.21 


1.17 


1.07 


±0.01 


±0.02 


xlO^ 


1.21- 


1.31 


1.26 


1.02 


±0.01 


±0.02 


xlO^ 


1.31- 


1.41 


1.36 


9.78 


±0.10 


±0.21 




1.41- 


1.52 


1.46 


9.38 


±0.10 


±0.20 




1.52- 


1.63 


1.57 


8.72 


±0.09 


±0.18 




1.63- 


1.76 


1.70 


8.59 


±0.09 


±0.18 




1.76- 


1.90 


1.83 


7.85 


±0.08 


±0.17 




1.90- 


2.04 


1.97 


7.41 


±0.07 


±0.16 




2.04- 


2.20 


2.12 


7.03 


±0.07 


±0.15 




2.20- 


2.37 


2.28 


6.38 


±0.06 


±0.13 




2.37- 


2.55 


2.46 


6.01 


±0.06 


±0.13 




2.55- 


2.75 


2.65 


5.45 


±0.05 


±0.13 




2.75- 


2.96 


2.86 


5.02 


±0.05 


±0.12 




2.96- 


3.19 


3.08 


4.62 


±0.05 


±0.11 




3.19- 


3.44 


3.32 


4.17 


±0.04 


±0.10 




3.44- 


3.71 


3.57 


3.79 


±0.04 


±0.09 




3.71- 


3.99 


3.85 


3.37 


±0.04 


±0.08 




3.99- 


4.30 


4.14 


3.02 


±0.03 


±0.07 




4.30- 


4.63 


4.47 


2.74 


±0.03 


±0.06 




4.63- 


4.99 


4.81 


2.41 


±0.03 


±0.06 




4.99- 


5.38 


5.18 


2.11 


±0.02 


±0.05 




5.38- 


5.79 


5.58 


1.87 


±0.02 


±0.04 





p 

(GcV/c) 



Flux±AFsta ± AFs, 



(m 



-i(GcV/c)-i) 



0.598 1.25 ±0.02 ±0.03 xlO^ 

0.644 1.22 ±0.02 ±0.03 xlO^ 

0.694 1.19 ±0.02 ±0.03 xlO^ 

0.748 1.15 ±0.02 ±0.02 xlO^ 

0.806 1.13 ±0.01 ±0.02 xlO^ 

0.868 1.08 ±0.01 ±0.02 xlO^ 

0.936 1.06 ±0.01 ±0.02 xlO^ 

1.01 1.00 ±0.01 ±0.02 xlQi 

1.08 0.96 ±0.01 ±0.02 xlO^ 

1.17 0.91 ±0.01 ±0.02 xlQi 
1.26 0.88 ±0.01 ±0.02 xlO^ 
1.36 8.53 ±0.10 ±0.18 

1.46 8.01 ±0.09 ±0.17 

1.57 7.53 ±0.08 ±0.16 

1.70 7.22 ±0.08 ±0.15 

1.83 6.72 ±0.07 ±0.14 

1.97 6.28 ±0.07 ±0.13 

2.12 5.71 ±0.06 ±0.12 

2.29 5.36 ±0.06 ±0.11 

2.46 4.92 ±0.05 ±0.10 
2.65 4.62 ±0.05 ±0.10 
2.86 4.09 ±0.05 ±0.09 
3.08 3.69 ±0.04 ±0.08 
3.31 3.31 ±0.04 ±0.07 

3.57 3.06 ±0.04 ±0.06 
3.85 2.70 ±0.03 ±0.06 
4.14 2.37 ±0.03 ±0.05 

4.47 2.11 ±0.03 ±0.04 
4.81 1.87 ±0.02 ±0.04 

5.18 1.66 ±0.02 ±0.04 

5.58 1.46 ±0.02 ±0.03 
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Table 3 

Atmospheric muon flux at sea level (continued). 



Momentum 

range 
(GcV/c) 



P Flux±AF,ta ± AFays P 

(GcV/c) (m-2sr-is-i{GeV/c)-i) (GcV/c) 



M 

FluxiAFsta ± AFsys 
(m-2sr-ls-l(GeV/c)-l) 



5.79- 
6.24- 
6.73- 
7.25- 
7.81- 
8.41- 
9.06- 
9.76- 
10.5- 
11.3- 
12.2- 
13.2- 
14.2- 
15.3- 
16.4- 
17.7- 
19.1- 
20.6- 
23.9- 
27.7- 
32.1- 
37.3- 
43.3- 
50.2- 
58.3- 
67.7- 
78.5- 
91.1- 
106.- 
132.- 
165.- 
207.- 
258.- 
323.- 



6.24 
6.73 
7.25 
7.81 
8.41 
9.06 
9.76 
10.5 
11.3 
12.2 
13.2 
14.2 
15.3 
16.4 
17.7 
19.1 
20.6 
23.9 
27.7 
32.1 
37.3 
43.3 
50.2 
58.3 
67.7 
78.5 
91.1 
106. 
132. 
165. 
207. 
258. 
323. 
404. 



6.01 
6.48 
6.98 
7.52 
8.10 
8.72 
9.40 
10.1 
10.9 
11.8 
12.7 
13.6 
14.7 
15.8 
17.1 
18.4 
19.8 
22.1 
25.6 
29.7 
34.5 
40.1 
46.5 
54.0 
62.7 
72.8 
84.4 
98.0 
118. 
147. 
184. 
229. 
286. 
359. 



1.63 ±0.02 ±0.04 
1.44 ±0.02 ±0.03 

1.21 ±0.02 ±0.03 

1.06 ±0.01 ±0.02 

9.11 ±0.13 ±0.21 xlO-i 

8.07 ±0.11 ±0.19 xlO-i 

7.06 ±0.10 ±0.16 xlO-i 
5.67 ±0.09 ±0.13 xlO-i 
4.90 ±0.08 ±0.11 xlO-i 
4.38 ±0.07 ±0.10 xlO-i 
3.58 ±0.06 ±0.08 xlO~i 

3.07 ±0.06 ±0.07 xlQ-i 
2.49 ±0.05 ±0.06 xlQ-i 
2.13 ±0.04 ±0.05 xlO-^ 
1.76 ±0.04 ±0.04 xlO-i 
1.49 ±0.03 ±0.03 xlO-i 

1.22 ±0.03 ±0.03 xlO-i 

9.57 ±0.09 ±0.22 xlO^^ 
6.69 ±0.07 ±0.16 xlO-2 
4.52 ±0.06 ±0.10 xlO-2 
2.93 ±0.04 ±0.07 xlO^^ 
2.01 ±0.03 ±0.05 xlO-2 
1.31 ±0.02 ±0.03 xlO-2 
8.76 ±0.18 ±0.22 xlO'^ 
5.72 ±0.14 ±0.15 xlO-'^ 
4.01 ±0.11 ±0.11 xlO-'^ 
2.55 ±0.08 ±0.08 xlO"'^ 

1.55 ±0.06 ±0.05 xlO-3 
9.43 ±0.43 ±0.40 xlQ-^ 
4.26 ±0.26 ±0.22 xlO""' 
2.51 ±0.18 ±0.15 xlO""' 
1.18 ±0.11 ±0.09 xlO-* 

5.58 ±0.68 ±0.52 xlO'^ 

2.56 ±0.41 ±0.30 xlO-5 



6.01 
6.48 
6.98 
7.52 
8.10 
8.72 
9.40 
10.1 
10.9 
11.8 
12.7 
13.7 
14.7 
15.8 
17.0 
18.4 
19.8 
22.1 
25.6 
29.8 
34.5 
40.1 
46.5 
54.0 
62.7 
72.8 
84.5 
97.7 
117. 
146. 
185. 
228. 
286. 
355. 



1.24 ±0.02 ±0.03 
1.13 ±0.02 ±0.02 
0.96 ±0.01 ±0.02 
0.83 ±0.01 ±0.02 
6.93 ±0.11 ±0.15 xlO-i 
6.08 ±0.10 ±0.13 xlO-i 
5.15 ±0.09 ±0.11 xlO-i 

4.54 ±0.08 ±0.10 xlO-i 
3.85 ±0.07 ±0.08 xlO"! 
3.24 ±0.06 ±0.07 xlO'^ 
2.74 ±0.06 ±0.06 xlO^^ 
2.31 ±0.05 ±0.05 xlO-i 
1.97 ±0.04 ±0.04 xlO-i 
1.68 ±0.04 ±0.04 xlO-i 
1.40 ±0.03 ±0.03 xlO-i 
1.15 ±0.03 ±0.02 xlO-i 
0.96 ±0.03 ±0.02 xlO-i 
7.36 ±0.08 ±0.15 xlO-2 
5.04 ±0.06 ±0.11 xlO-2 
3.48 ±0.05 ±0.07 xlO'^ 
2.26 ±0.04 ±0.05 xlO-^ 

1.55 ±0.03 ±0.03 xlO-2 
1.04 ±0.02 ±0.02 xlO-2 
6.91 ±0.16 ±0.16 xlO-3 
4.54 ±0.12 ±0.11 xlO-3 
2.88 ±0.09 ±0.08 xlO'^ 
1.93 ±0.07 ±0.06 xlO-3 
1.17 ±0.05 ±0.04 xlO-3 
6.95 ±0.37 ±0.28 xlO""* 
3.65 ±0.24 ±0.18 xlO"* 
1.47 ±0.14 ±0.09 xlO"* 
0.75 ±0.09 ±0.06 xlO"'' 
4.13 ±0.58 ±0.39 xlO'^ 
2.45 ±0.40 ±0.28 xlO'^ 
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